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In the chiral magnetic effect, there is a competition between a strong magnetic field, which 
tends to project positively charged particles to have spin aligned along the magnetic field, and a 
chirality imbalance which may be produced locally by a topologically nontrivial gauge field such 
as an instanton. We study the properties of the Euclidean Dirac equation for a light fermion in 
the presence of both a constant abelian magnetic field and an SU(2) instanton. In particular, we 
analyze the zero modes analytically in various limits, both on R* and on the four-torus, in order to 
compare with recent lattice QCD results, and study the implications for the electric dipole moment. 
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I. INTRODUCTION 

Since quarks carry both electric and color charge they couple to both electromagnetic and gluonic gauge fields. A 
magnetic field introduces a Landau level structure to the fermion spectrum, in which the zero modes have definite 
spin, aligned along the magnetic field [TH5]. In a gluonic field with nontrivial topological charge the fermion spectrum 
also has zero modes, with chiralities determined locally by the local topological charge of the gauge field [5l413j. The 
associated instanton transitions in the 0- vacuum can lead to a fluctuating electric dipole moment of the neutron jl4j . 
In this paper we investigate what happens when a quark experiences both a strong magnetic field and a topologically 
nontrivial gluonic field, such as an instanton. For a single instanton the fermion spectral problem has a conformal 
symmetry [151 1161 . and the zero modes are localized on the instanton, falling off as a power law with Euclidean 
distance. The conformal symmetry is broken by the introduction of a magnetic field, and now the zero modes develop 
an asymmetry, falling off in Gaussian form in the plane transverse to the B field, but as a power law in the other two 
directions. This basic asymmetry is the key to the phenomenon of magnetic catalysis |17j and the chiral magnetic 
effect [TBHH], as sketched in Fig. 1. In this paper we discuss some features of the spectral problem for fermions 




FIG. 1: A sketch of the topological charge density, q oc trJ^fj,,^jFi^^, for a single instanton [red], and the density of the quark 
zero mode [grey]. On the left, there is a single instanton, and both densities fall off as power laws, with q falling off faster. On 
the right, with the introduction of a magnetic field, the topological charge density is unchanged but the zero mode density is 
distorted into an asymmetric shape, localized along the direction of the strong magnetic field. 
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in the combined background field of a magnetic field and an instanton. We are motivated by situations in which 
quarks experience both types of fields, such as in dense astrophysical objects such as neutron stars and magnetars, 
and in heavy ion collisions such as those at RHIC and at CERN [5DJ [231 HI] • We are also motivated by recent lattice 
QCD analyses |25H29| . which provide important numerical information about the Dirac spectrum in both QCD and 
magnetic field backgrounds. Analytically, while the effect of each individual background is very well known, their 
combined effect turns out to be quite intricate. 



II. GENERAL FORMALISM: DIRAC SPECTRUM 



We briefly review the well-known properties |13j of the Dirac equation for an instanton background or a constant 
magnetic field, since our goal is to discuss what happens when we combine the two background fields. We work in 
Euclidean four-dimensional spacetime, with the following conventions. We follow the notation of [TD], and express the 
4x4 Dirac matrices, 7^, for = 1, 2, 3, 4, in terms of the 2x2 matrices — (1, —ia) and = (1, ia) = aj,, [here 
a are the usual 2x2 Pauli matrices]: 

^ (5, "o") ' ^5 ^ (0 -1) 

Thus, the Euclidean Dirac operator can be expressed as 



a,,V^\ f D 



a,,V^ j -\-D^ Q 



(2) 



where the covariant derivative, — — is written with a hermitean gauge field, A^, and is the Euclidean 
time coordinate. We write the gauge field Af^ as a sum of a non-abelian part, A^i, and an abelian part, a^: 

A^ = A^, + (3) 

with the respective coupling constants absorbed into the gauge fields. The Dirac operator is anti-hermitean, so we 
write (with A real) 

i'pxPx = \^x (4) 

Since {75,?^} = 0, we can take A in Q to be non- negative, with the negative eigenvalue solutions simply given by 
ip-x = "fsipx- This means that we can effectively discuss the zero modes (A — 0) separately, and for the nonzero modes 
(A 7^ 0) we consider the squared operator: 

m'i^x^C"^' n?nV^A = AVA (5) 



\^ D'fD 

The positive chirality sector, x = +1; is described by the operator DD\ while the negative chirality sector, x — ~li 
is described by the operator D'^D. We can write these operators as 

X = +l: DD'^ = -Vl-T^^G^,, (6) 
X = -l: D^B = -Vl-F^,a^, (7) 

We have used \D^^,'Di,] — —iJ-f^u, where J^^^ is the field strength associated with the gauge field Af^, and the spin 
matrices (j^jy and a^^ are defined as 

^ni' = ^ (a^aiy - ai,af_i) (8) 
<^fii^ = ^ (a^az. - a^a^) (9) 

In (6[7) we have used the properties PJ]: Ouai, — Sf^^, + 2i(T^j^, and a^ajy = (5^,^ -I- 2ia^jj. 

For non-zero modes [i.e., solutions to ([5| with A 7^ 0], the operators DD'^ and D have identical spectra, for 
any background field. This is simply because we have an invertible map: suppose the 2-component spinor v satisfies 

Dv = X^v. Then u = Dv is clearly an eigenfunction of the other operator, DD\ with precisely the same eigenvalue: 
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DD^u = DD^ Dv = X^u. Similarly, if u satisfies DD^u = X'^u, then v = D^u is an eigenstate of D^D with the same 
eigenvalue. Thus, when A 7^ 0, we can write the 4-component spinor solution in the form 

= f_i'55tuj ^here DD^u^^X'ux (10) 



A 



or in the form 



This is true for any background field: non-abelian, abelian, or both 



^>^=i^y^) ^he^e D'^Dvx = X^vx (11) 



A. Magnetic field background 

For a constant (abelian) magnetic field, of strength B, pointing in the X3 direction, we have an abelian field strength 
/12 = B, and so we find 

X = +l: DD^ = -Vl-Bas (12) 
X = -1 : D^D ^ -Vl - Bas (13) 

where we have used the fact that CT12 = 1712 — ^cr^. Note that in this case the 2x2 operators DD^ and D^D, of the 
two chiral sectors, are the same, and therefore they have identical spectra. Due to the subtraction term, —Ba^, it is 
possible to have zero modes, and since DD^ = D^D these zero modes occur in each chiral sector. More explicitly, we 
can make a Bogomolnyi-style factorization and write 

-Vl-Bas - ~dl-dl-iViTi-D2)iV,±W2)±B-Ba3 (14) 
= - 9| - V^V± ±B-Ba3 (15) 

For zero modes, we take 93 = 94 = 0, and with i? > we choose the upper signs to ensure normalizable modes. For 
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FIG. 2: A sketch of the form of the Landau level spectrum of the squared Dirac operator {tPY , for a constant abelian 
magnetic field. The plot is the same for both positive and negative chirality; that is, for the operators DD^ and D. Note 
that for each chirality the zero mode has spin up (along the magnetic field), and furthermore note that the two-dimensional 
operators T>+T>- and are isospectral apart from the zero mode. 



example, in the symmetric gauge where the abelian gauge field 

B 



2 (-a;2,a;i,0,0) (16) 
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the zero modes can be expressed in terms of the normaHzable solutions to {Di + 0^2) u = 0: 



V'o = 9izi)e' 



or = 5(^1)6' 



-B\z,\^/2 



(17) 



Here 3(2:1) is a holomorphic function of the complex variable zi = {xi + 1x2) /V^- Both sets of zero modes have spin 
up, aligned along the B field; this is just the familiar lowest Landau level projection onto spin up states. Note also 
that the zero modes have the characteristic Gaussian factor in the {xi,X2) plane, transverse to the direction of the 
magnetic field. This factor is the origin of the distortion sketched in the right frame of Fig. 1. 

The number of zero modes per unit two-dimensional area [in the (xi, X2) plane] is given by the Landau degeneracy 
factor, the magnetic flux per unit area: B/{2t:). In fact, even for an inhomogeneous magnetic field B(xi, X2), pointing 
in the X3 direction, the number of zero modes [of each chirality] is determined by the integer part of the magnetic 
flux (this is the essence of the Aharonov-Casher theorem [Tj). For example, on a torus [5]: 



N+ = = 



1 

2^ 



d^xB 



(18) 



The higher Landau level states are the same for both spins, as {—'D_'D+ + B) and (— 1?+2?_ — B) have identical 
spectra, apart from the lowest level, which only has spin aligned along the magnetic field. The resulting spectrum is 
sketched in Fig. [2j 



B. Instanton background 

For an instanton field, A^, the (non-abelian) field strength i^^^ is self-dual [that is: F^j^i, — F^^, where the dual 
tensor is defined: F^^, = ^e-fiua/iFap]. Then the anti-self-duality property of a^jy [that is: a^^i, — —<J^u] implies: 

X = +1 : DD^ = -Vl (19) 
X = -l: D^D = -Vl-F^,a^, (20) 

Since —2?^ is a positive operator, this means that for an instanton background there can be no zero mode in the 
positive chirality sector. On the other hand, due to the subtraction term, —Ffj^^a^i^, in D^D, it is possible to have a 
zero eigenvalue solution in the negative chirality sector, and it has the form 

V'o , where Dv ^ (21) 

[For an anti-instanton, an anti-self-dual field with _F^j^ = —F^^, the zero mode lies in the positive chirality sector, 
because tr^i, is self-dual: a^n = a^u.] For a general non-abelian gauge field A^, which is neither self-dual nor anti-self- 
dual, the Atiyah-Singer index theorem [13 [H] states that the difference between the number of positive and negative 
chirality zero modes is given by the topological charge of the gauge field: 



N+-N^ = -^ J d'x F^,F^, (22) 

Here we have written Fi^^, = F^^T"", with generators normalized as tT{T°'T^) = ^(5°^. For gauge group SU{2), with 
fermions in the defining representation, we take generators T"^ = ^r" in terms of the Pauli matrices t, and we can 
write the single instanton gauge field |30j . centered at the origin, in the regular gauge as 

^ x^ + 

where p is the instanton scale parameter, and ry^^, is the self-dual 't Hooft tensor [SI [TU]. The topological charge 
density is 

1 ~ 192 

'^^''^ " 32^^ ^^"^^^^ " {x^ + p2)4 (24) 
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FIG. 3: A sketch of the form of the spectrum of the squared Dirac operator (i'P)^, for an instanton field. Note that the 
operators DD^ and D are isospectral [also with —'D^] except for a zero mode in the negative chirality sector. 

There is a single zero mode [6H9l I13j. also localized at the origin, with density: 



I'^ol' = t4^A^ (25) 

These densities both fall off as power laws, with scale set by p, but the topological charge density is more localized, 
as indicated in the left-hand frame of Fig. [l] The nonzero modes are given by (10) or (111, and we note that the 



spectra are identical in each chiral sector, apart from the zero modes, as is indicated in the sketch in Fig. [3j 

C. Combined instanton and magnetic field bacicground 

Physically, an instanton field projects the zero modes onto a definite chirality, while a constant magnetic field 
projects the zero modes onto definite spin, aligned along the direction of the magnetic field. When we combine the 
two background fields, both a non-abelian instanton field -F^j, and an abelian magnetic field /12 = S, there is a 
competition between the two projection mechanisms, and the outcome depends on their relative magnitude, as we 
show below. Technically speaking, the instanton zero mode has a specific ansatz form that unifies space-time and 
color indices, while the magnetic zero modes have a natural holomorphic structure, and these two different ansatz 
forms do not match one another. The competition between these two ansatz forms makes the combined problem 
nontrivial. For an instanton field, since the field falls off as a power law, all eigenmodes also fall off with power law 



behavior. On the other hand, once a constant magnetic field is introduced, for example in the gauge (16), all the 
eigenstates (even those in the higher Landau levels) have a Gaussian factor exp(— i?|zip/2) that localizes the modes 
near the axis of the magnetic field. This is the reason for the distorted density in the right-hand frame of Fig. 1. In 
the extreme strong magnetic field limit this leads to a dimensional reduction to motion along the magnetic field, with 
interesting physical consequences such as magnetic catalysis \TT and the chiral magnetic effect [T8', '2T| l22] . 



Concerning zero modes, we begin with a simple but important comment: in the index theorem (22), the magnetic 
field makes no contribution, since with the field strength decomposed into its non-abelian and abelian parts, J^^^ ~ 
F^_,y + /p,., we have 

tr (j-^, J"^.) = tr (i^p.F^.) + (dim) /^,/^, (26) 
= tr(F^,f'^,) (27) 

where dim is the dimension of the Lie algebra representation of the non-abelian gauge fields. The cross terms vanish 
since the Lie algebra generators T" are traceless, and the f^uffiv term vanishes since there is no abelian electric field. 
For example, if there is no n onabelian field, just an abelian magnetic field, then the topological charge clearly vanishes. 



and the index theorem (22) is consistent with the fact that DD^ = D for an abelian magnetic background (recall 
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(12 13 1), so that there are the same number of zero modes in each chiral sector. Now, with both background fields 



present, we find 

DD^ = -Vl-Ba^ (28) 
D^D ^ -Vl~F^,<j^,~Baz (29) 

Notice that the eigenvalues of DD^ are simply those of the scalar operator —2?^, with a spin term ±i3, as can be 
seen clearly in Figure [s] The fact that there is a subtraction term from the positive operator —2?^ in both chirality 
sectors tells us that it is possible to have zero modes for each chirality, but their number will depend on the relative 
magnitude of F and B. In the next Section we study a specific model where we can quantify this precisely. Another 
important implication is that we may also have some " near-zero- modes" , where the F and B subtractions do not 
exactly cancel the lowest eigenvalue of —2?^, but lower the eigenvalue of DD^ or D^^ D to near zero. 



III. LARGE INSTANTON LIMIT 

In the very strong magnetic field limit, where the magnetic length, 1/y/B, is small compared to the instanton size p, 
we expect a significant distortion of instanton modes and currents. In this limit we can make a simple approximation 
that reduces the problem to a completely soluble system. 



A. Covariantly constant SU{2) instanton and constant abelian magnetic field 

In the large instanton limit, we expand the instanton gauge field as: 

« '^'n>'^ + • ■ • (30) 

To leading order in such a derivative expansion, the non-abelian gauge configuration A^(x) is self-dual and has 
covariantly constant field strength: F^^, = rj'^^. In this limit we can make an SU{2) "color" rotation, along with 
a choice of Lorentz frame, to make the instanton field diagonal in the color space (we choose the direction), so 
that the field is self-dual, covariantly constant and quasi-abelian. Defining the instanton scale F = the combined 
gauge field, including also the abelian magnetic field as in ([S]), can be written as: 

F B 

Af, = - — {-X2,Xi,-X4,X3)t^ + —(-X2,Xi, 0,0)12x2 (31) 

This gauge field is fully diagonal and moreover is linear in x^, so the problem is analytically soluble (this is the basic 
premise of the derivative expansion). The only nonzero entries of the field strength tensor are 



J-,, = -Ft^ + BI 



B- F 
B + F 



-^34 = -Ft' = ( (f ^p) (32) 

In the absence of the magnetic field the field strength is self-dual, J- 12 = ^34, but a nonzero magnetic field breaks this 
symmetry. The topological charge density is (recall the normalization of the generators) 



-pa -pa _ _\ /_ _ /oo\ 

n.v' uv or,_2 0_2 ^^^1 



32^2 327^2 27r2 

To study the Dirac spectrum we consider the 2x2 operators DD"^ and D'^D in (6[7). Notice first that 



^iiv CTpi/ = (-7^12 - -^34) cr3 (34) 
^iiv^^iiv = (-7^12 + -^34) cr3 (35) 

It is convenient to factor the 4-dimensional Euclidean space and consider separately the (a:i, X2) plane and the (^3, 0:4) 
plane, as sketched in Figure |4j Then in the (xx^x^) plane we have a (relativistic) Landau level problem with effective 
field strength (B — F) in the = -1-1 sector, and with effective field strength {B + F) in the — —1 sector. In the 
(0:3, X4) plane we also have a (relativistic) Landau level problem, now with effective field strength —F in the — +1 
sector, and with effective field strength F in the = —1 sector. In the {xi,X2) plane the sign of the effective field 
strength depends on which of i? or F is larger, so we consider separately the cases B > F or B < F. 






FIG. 4: Sketch of the effective magnetic field strengths perpendicular to the {xi,X2) and {x3,X4,) planes and the associated 
zero mode densities which are represented by the colored disks. The colors denote the isospin dependence: blue (upper) stands 



for T = 1, and red (lower) stands for r = —1. All the zero modes are spin up {a = 1) for B>F. See eqn. (61 1 for the case of 
4-torus where the magnetic fluxes are quantized and coincide with the number zero modes. 



1. Strong magnetic field limit: B > F 

When B > F, both {B — F) and {B + F) are positive. Thus, each color component of J^i2 is associated with a 
positive "magnetic" field. On the other hand, for 7^34, the — +1 sector has a negative field strength, while the 
— ~l sector has a positive field strength. 

We first consider the = +1 case. Then J^i2 — [B — F), 7^34 = —F, Ff^^a^y = Bay,, and F^„(jf^y — [B — 2F)a3. 
With a positive field strength the normalizable zero state is given by (2?i + il?2) u — 0. But since J-34 is negative, 
we factorize the corresponding covariant derivatives in the opposite order, in order to obtain a normalizable state 
annihilated by (I?3 — iV^). Thus, we have: 



X = +1 : DD'f = 
X = -l: D^D = 



(2?i - iV2) {Vi + iV2) + F12 - {V3 + iVi) [Vs - iVi) - Fsi - Bay 

(2?i - W2) (Pi + iV2) - (2?3 + i'Di) (Vy - Wi) + B-Bcji (36) 

(2?i - W2) {Vi + iV2) + F12 - {Vy + iVa) (I?3 - i'Di) - -^34 -{B- 2F)ay 

(2?i - W2) (Pi + iV2) - {V-3 + iVi) (P3 - iVi) +B-{B~ 2F)ay (37) 



This shows that there is a zero mode, when the spin term Bay cancels the B term from the Bogomolnyi factorization 
of the covariant derivative term. This occurs in the positive chirality sector, x = +lj ^J^d with spin up: ay = +1. 



Now consider the = —1 case. Then F12 
field strengths are positive, so we write 



X = +1 : DD'^ = 
X = -1 : D'^D = 



[Vi - iV2) {Vi + tV2) 
{Vi - iV2) {Vi + W2) 
{Vi - iV2) {Vi + W2) 
{Vi - iV2) {Vi + W2) 



{B + F), J-34 = F, J-^,a^, = Bay, and -F^.tr^. = {B + 2F)ay. All 

Fl2 - {Vy - iVi) [Vy + iVi) + J-34 - Bay 

{Vy - iVi) {Vy + iVi) + {B + 2F) - Bay (38) 

Fi2 - {Vy - 1V4) {Vy + iVi) + J-34 - {B + 2F)ay 

{Vy - iVi) {Vy + iVi) + {B + 2F) -{B + 2F)ay (39) 



This shows that there is a zero mode, but now in the opposite chirality sector, x = ^Ij £^nd also with spin up: 
oy = +1. 

To summarize: when B > F, the Ty ~ +1 color sector has spin up zero modes with positive chirality, while the 
T3 = — 1 color sector has spin up zero modes with negative chirality. We can count the number of zero modes in each 
chirality sector by simply taking the product of the Landau degeneracy factors for the {xi,X2) and (x3,a;4) planes, 
with the corresponding effective magnetic field strengths. Therefore, the corresponding Landau degeneracy factors 
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FIG. 5: Spectrum of the sq uare d D irac ope rator (ip)"^ , for both a strong magnetic field and instanton background, with B ^ F, 
as derived from equations ([36] [37) |38] 39|. Note that for each "color", T3 = ±1, the operators DD^ and D are isospectral 
except for a zero mode. For T3 = +1, this zero mode is in the positive chirality sector, while for T3 = —1, this zero mode is in 
the negative chirality sector. The spin-projection steps relating the spectra of DD\ D and DjJ, steps of ±_B, ±(5 — 2F) 
and ±(-B + 2F), are indicated in the figures. 



give the zero- mode number densities (i.e, the number per unit volume): 

(B-F) F 



X = +1 : 
X = —1 : ri- 
The index (density) is the difference, 



2tt 2-k 

{B + F) F_ 
2tt 2-k 



(t3 +1 , erg = -hi) 
(1-3 = -1 , ff3 = +1) 



(40) 

(41) 



n I — n_ = — 



27r2 



(42) 



in agreement virith the general index theorem (22 1, in view of (33). We also note that the total number density of zero 
modes 



BF 
2^ 



(43) 



is linearly proportional to the magnetic field strength B. This is in agreement with numerical lattice gauge theory 
results [Ml. 



2. Weak magnetic field limit: B < F 



Even though this limit is outside the derivative expansion regime that originally motivated the quasi-abelian co- 
variantly constant gauge field ansatz (30), it is still instructive to see how the Dirac spectrum changes for such a 
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field when B < F. The difi'erence is that now the = +1 effective "magnetic" field {B — F) is negative, so the 
factorization should be done in the opposite order also for the (xi, X2) plane. 

Consider the = +1 case first. Then J^i2 = (B — F), = —F, F^^d^^, — -60-3, and Ff^^a^i, = [B — 2F)<7^. Since 
F12 has the opposite sign from before, we factorize the corresponding covariant derivatives in the opposite order, in 
order to obtain a normalizable state annihilated by {T>i — il?2)- Thus, we have: 

X = +l: DD^ = - {Vi + W2) (Vi ~ iV2) ~ F12 - {V3 + tV4) {V3 - iVi) - F34 - Ba3 

= - (2?i + iV2) (Pi - iV2) - (P3 + il^i) (2?3 - i2?4) -{B- 2F) - Ba^ (44) 

X = -l: D^D = - (Vi + W2) {Vi - tV2) - F12 - {V3 + tVi) {V3 ~ W4) ~ F34 - {B - 2F)a3 

= - (2?i + iV2) (Pi - tV2) - (P3 + iVi) {V3 - iV^) -{B- 2F) - {B - 2F)a3 (45) 

This shows that there is a zero mode, in the negative chirality sector, x — — 1, and with spin down: (T3 = —1. 

Now consider the = — 1 case. This is exactly as before, with F12 = {B + F), F34 = F, Ff^^a^i, = Ba^, and 
^fj.uO'^v — {B + 2F)(j3. Thus, we can write DD^ and D'^ D exactly as in (38 39), and we see that, as before, the zero 
modes are in the negative chirality sector, with spin up. 

To summarize: when B < F, the T3 = +1 color sector has spin up zero modes with negative chirality, while the 
T3 = — 1 color sector has spin down zero modes also with negative chirality. Counting the corresponding Landau 
degeneracy factors we obtain: 

X - +1 : n+ = (46) 

" (B+F) F / -, „ I 1^ 

2¥ ' (^3 = -1 , <y3 = +1) 

(t3 = +1 , (73 = -1) 



27T 2tt 

-B+F) F 



2tt 2tt 

The total number density of negative chirality zero modes is (note that the B dependence cancels) 

n- = - — 
27r 

As before, the index (density) is given by the difference, 

n,-n_=-|J (49) 



(48) 



in agreement with the general index theorem ( 22 ) . In this case the total number density of zero modes 

= (50) 
which is independent of the magnetic field strength B, and equal to (minus) the index. 



B. Physical picture 

These results lead to the following simple physical picture. The instanton tries to generate a chirality imbalance but 
is neutral to the spin, whereas the magnetic field tries to generate a spin imbalance but does not affect the chirality. 
Depending on which is stronger, the zero modes have either a definite spin with a chirality imbalance (B > F), or a 
definite chirality with a spin imbalance {F > B). Also we see that in the former case, the total number of zero modes 
scales with B and is not equal to the index, while in the latter case it is independent of the B field, and is equal to 
the magnitude of the index (even though the field is not self-dual). 

More explicitly, for the B > F case, consider starting with just a strong magnetic field B, later turning on a weak 
instanton field. Without the instanton field, the zero modes and their degeneracy are given by the Aharonov-Casher 



theorem ( 18 ), so that the zero mode density is the Landau degeneracy factor B/(2t:) for each chirality sector. All the 
zero modes are spin up, as is familiar for the lowest Landau level (see Fig. [2]). There is an equal number of positive 
and negative chirality zero-modes, which is consistent with the index theorem, since the topological charge vanishes 
for a constant B field. Now consider turning on an instanton field F, with B > F > 0. We see from ( 40|41 ) that 

the effect of the instanton is to flip some of the chiralities: (^) positive chirality modes become negative chirality 
modes, leading to a chirality imbalance of in agreement with the index theorem (22). On the other hand, the 
total number of zero modes, grows linearly with the magnetic fleld when F is nonzero. 
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f2 

On the other hand, when F > B we have the following physical picture. With just the instanton, there are 

zero modes, all with negative chirality. Of these, (^)^ have spin up, and (^)^ have spin down. Now turn on a 
magnetic field B, with < B < F. The eff ect of the magnetic field is to flip some of the down spins to up spins, 
without affecting the chirality. From (46 47 1 we see that ^^^2 zero modes have their spin flipped, leading to a spin 
imbalance, without creating a chirality imbalance. Thus, the index is still equal (in magnitude) to the total number 
of zero modes. 



C. Creation and Annihilation Operator Formalism on the Four- Torus 

To count degeneracies it is convenient to introduce torus boundary conditions, ~ + L^, and so to study this 
configuration on a 4-torus. The nontrivial homotopy group of the torus leads to a well defined topological structure 
for a constant instanton field [32H35j . Furthermore, it allows a direct comparison with some recent lattice results 
^31171125], in particular the counting of zero modes. We use torus boundary conditions, namely that the gauge field 



(31 1 is periodic up to a gauge transformation: 

A^{x^ + L^) = n-^{x){A^,{x^) - id^,)Q.^{x) (51) 

where v denotes the shifted coordinate and is the associated cocyle. A general treatment of non-abelian gauge 
fields on the torus can be found in [32H35] . 

As in the standard analysis of the constant magnetic field problem, it is useful to work with complex coordinates. 
The 4-torus can be imagined as two orthogonal 2d planes, with appropriate periodicity conditions, parameterized by 
two complex coordinates: 

Xi+ix2 xz + ixi 

In these coordinates, the 2x2 Dirac operator D = a^V^, defined in ([2j), is 

d, - ^-zA (82 + f Z2 9l - ^^-1' 




D = -iV2{ \(»i+-iV2\ |«)D_ (53) 

-92 - f ^2 / \dl + -82 + f Z2 

Here t± denote the color projection matrices. Following the standard treatment for a constant magnetic field, we 
define the ladder operators: 

ai ^ —iV2 ( Bi + ^-—^zi] , for D+ 



- B I F 
ai = —i\f2 I d\ H — zi ) , for 



02 = -«%/2 ^(?2 + -^^2^ , for D+ 

~a2 = -i^{d2 + ^z^ , forD_ (54) 
These satisfy the commutation relations 

[ai,4] =2(B-F) 
[ai,al] =2(B + F) 
[a2,4] = 2F 

[a2,4] = 2F (55) 

It is now explicit that, for each color, we have a set of two independent Landau level problems and the Landau levels 
are governed independently by the annihilation-creation operators (ai,a|) for T3 = -1-1, and (oi,aJ) for ''"3 = We 
denote the associated number operators as: 

4ai ==2(B-F)iVi , a\a2^2FN2 

ajoi = 2(B + F)iVi , 4a2 = 2F7V2 (56) 
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One should keep in mind that within our large instanton approximation, in which B ^ F, the Landau levels for the 
fields perpendicular to the {xi,X2) plane have greater degeneracy than the ones for the fields perpendicular to the 
{x3,X4) plane. Expressed in terms of these ladder operators, the Dirac operator D = a^V^, defined in ([2]), is given 
by: 

«n^D +p2 «1'\^„_ (57) 
\ai -a2j \ai -a\J 

To find the Dirac eigenvalues we consider {ifiY , as in ([5|, which means we require: 

DL>^ = ( 2(B - F)Ni + 2FN2 + f [! „ „ ) ) «) D+ + (2(5 + F)Ni + 2FN2 + ' ° 



2S;; ^ "+ ^ 1^"^^ ..2^ 2{B + F) 

D^D = (^2{B~F)N^+2FN2+{^^ 2(5 - F)) ) ® "+ + (^(^ ^ + 2^^2 + (^0 2(i3 + 2F)))®"- 

(58) 

Thus, {iPY is fully diagonalized and we have a complete description of the entire spectrum in terms of elementary 



harmonic oscillator number operators. Notice that, in addition to the Landau level operators (54), we can define 
another set of ladder operators which commute with all the and a^, and therefore with the Dirac operator ij): 

hi = -^\/2 - ^^^z^ , for D+ 

61 = -iV2 (^Bi - ^^^-^^i) , for 

62 = -iV2 ^82 - -^^2^ , for D+ 

62 = -iV2 (^32 - ^22^ , forD_ (59) 

These operators satisfy the same commutation relation as the Landau level ladder operators: 

{bi,b{]=2iB^F) 
[bi,b\]=2iB + F) 
[62,4] =2F 

[b2,bl]^2F (60) 

These operators generate magnetic translations and so characterize the degeneracies of the Landau levels. With torus 
boundary conditions, the degeneracy is finite and given by the net flux quanta through the period parallelogram of 
the corresponding 2-torus. In other words, each Landau level is a finite dimensional representation of the magnetic 
translation group [5J [51 135] with the dimension being equal to the flux. For simplicity of notation we take all four 
periods to be equal, denoted by L. Then for each of the two 2-tori [and for each of the two colors] we have the flux 
quantization conditions: 

{B - F)L^ =^ 2n{N - M) , for D+ 
[B + F)L'^ ^2t:{N + M) , for D_ 
FL^^2ttM , for D+ 

FL^^2ttM , for D_ (61) 

Here N and M are nonzero integers and each Landau level has degeneracy factor 

N±=M\NtM\ , for T3 = ±l (62) 

since we simply multiply the degeneracy factors for the two independent 2-tori for the (xa, X4) and (xi, X2) planes. The 
explicit states can be generated by acting on the ground state with the various creation operators, and the associated 
wavefunctions can be written in terms of elliptic functions, as is familiar PJH]. We will concentrate instead on the 
full eigenvalue spectrum of {ip)"^ . 



12 



D. Zero modes 

As described above, we treat separately the cases where B > F and B < F. The torus boundary conditions allow 
us to express the degeneracies exactly since the relevant fluxes are integers. 



1. B> F 



As explained above, with just the magnetic field present, there are only spin-up zero modes with both chiralities. 
Let us denote the degeneracy as 2iV, where BL^ /2tt = N. The factor 2 is due to color. Now consider turning on an 



instanton field F, so that B > F > 0, with F = 2ttM/L'^ as in 

N2 = Ni^ N2 = 0. It is clear from 



looking at (58), with A^i 

zero modes. Their spin and chirality can be read off directly. Furthermore, their degeneracies are fixed by (61 1 



. The zero modes can be easily constructed by 
that there only two sectors where one can have 

As a 



result, we see that there are (N — M)M spin up, positive chirality zero modes for the D+ color sector, and {N + M)M 
spin up, negative chirality zero modes for the D_ color sector. Here M is the instanton flux and N is the magnetic 
flux, from (61). The index of the Dirac operator is given by the difference: 



index (|!)) = N+- 



{N - M)M -{N + M)M 

F2 i4 



27r2 



(63) 



in agreement with (22) and (33). The effect of the instanton is to flip the chirality of Af^ fermion zero-modes, resulting 
in a chirality difference of 2AP. Also, the total number of zero modes is: 



total number of zero modes = 



[N + M)M + {N - M)M 
2NM 
BFL^ 
27r2 



(64) 



in agreement with ( 40|41 ). We see that the total number of zero modes is linearly proportional to the magnetic field 
strength. This agrees with the lattice results [28]. The functional forms of these zero modes can again be constructed 
in the same manner as above. 



2. F > B 



It is a straightforward exercise to construct the analogue of (58 ) for F > from which the zero modes can easily 
be constructed. First, consider just the instanton field, so F is nonzero, but B = 0. In this case, there are 2M^ zero 
modes. All of them are positive chirality, in agreement with the vanishing theorem [9] that states that for a self-dual 
field the index is equal to the total number of zero modes, since all zero-modes have the same chirality. Also, as 
expected, there is no preference in spin and there are equal number of spin up and down zero modes. 

When we turn on a weak magnetic field, so that B = 2ttN/L'^ is nonzero, the effect of the magnetic field is to flip 
some of the down spins to up spins. Now there are {N + M)M spin up zero modes, and {M — N)M spin down zero 
modes. The zero modes are still positive chirality, since the magnetic field does not flip chirality. Therefore the index 
is still equal to the total number of zero modes even though the gauge field is not self-dual anymore with the magnetic 
field. 

index(l^) = N+-N^ ^ - 2M'^ 
= -2M2 
F2 i4 



27r2 



(65) 



In this case, the total number of zero modes is also equal to the index: 



total number of zero modes 



(M - N)M + (M + N)M 
2M^ 



(66) 
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As a result we explicitly get the same simple physical picture as described above. Depending on whether the magnetic 
field or the instanton is stronger, the zero modes have either a definite spin with a chirality imbalance {B > _F) or a 
definite chirality with a spin imbalance {F > B). The advantage of the torus construction is two-fold: first, the fluxes 
are integers and the wavefunctions have simple expressions in terms of elliptic functions; and second, the counting of 
zero modes agrees with the numerical lattice QCD results. 



E. Landau levels 



From the quantum mechanical supersymmetry relations (10 11) we can construct the full spinor solution from 
l^fl), a two-component spinor satisfying: 

(67) 

Since DD^ commutes with cr^, from (29), we can choose to be either spin up or spin down. In the T3 = +1 
sector, we find 



ni,n2, 




A„,,„, = ^2{B-F)ni+2F 



n2 





ni - 1,^2 - 1; 



A„,.„, - ^/2{B - F) m + 2F 



n2 



(68) 



For the spin up states, ni, 71,2 G {0, 1, 2, 3...}, excluding the case where both ni = 712 = 0, while for the spin down case 
'ni,n2 € {1,2,3...} label the Landau levels corresponding to independent magnetic fields in the {xi,X2) and {x3,X4) 
planes, and |ni,n2) are normalized harmonic oscillator eigenstates. One should keep in mind that (for B > F) each 
level has a degeneracy M{N — M), with M and {N — M) being the fluxes in the (a;i,X2) and (0:3, 0:4) planes. The 
full spinors are: 



IV't) 



/ |ni,n2) 
1 

;7f -l^2Fn2\nun2-l) 

\-i^2{B-F)n,\n,^l,n2)J 

( \ 

\n\ - l,n2 - 1) 

y2{B-F)n^\n^ - 1,712 



(69) 



V -i^/2F7^^|nl- 1,712) / 

The overall factor fixes the normalization of the spinor to 1. In the T3 = — 1 sector, there is a similar construction, 
recalling that the field strength is (B -I- F) in the {x\^X2) plane. 



IV. MATRIX ELEMENTS AND DIPOLE MOMENTS 



In this Section we consider certain matrix elements involving quark bilinears, such as have been computed on the 
lattice. For these purposes, it is convenient to introduce a small quark mass 771, so that the propagator of the Dirac 
operator ^ + 711 is given by: 



1 / -Ti ^1 n 

1 nt ^ik 



(70) 



Note that DD^ and D have identical spectra, except for possible zero modes, so they can be viewed as square 
operators (matrices) of different dimension, as is clear when they are diagonalized in their respective eigenspaces. The 
zero mode contribution to the propagator can be separated by writing it in one of two ways, depending on which 
chirality supports zero modes 



-1 D 



1 nt m 

m^+D^D^ rrfi+DID 



(71) 
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An important set of quark bilinears involve the spin tensor Emi/ 



= ^[7.,7.]=2(Y f ) (72) 



This representation makes clear the natural decomposition of E^j^ into its self-dual part ((T^t,y) and its anti-self-dual 
part {o^v). The bilinears are 



(^S,,V)=tr(^S,.^^j (73) 
For applications to the chiral magnetic effect, we are interested in the magnetic and electric dipole moments: 

= \e^M^^Jk^) (74) 

af = (V;E,4V') (75) 
With a strong magnetic field in the X3 direction, we concentrate on and ct|', which require the spin tensors: 



El 



2 



fig 
fT3 



Thus, 

- / Tfl^ \ ( TfL^ \ 

mmi^i,) = tr..2 (^-3^,_^^^,j+tr2x2 (^'-3^,^^,^j (77) 

m(V;E34^) = -tr,x2 (^^3 ^,^^^, j + tr.x2 (^^3 ^,^^,^ j (78) 

The dominant contribution to the trace over the spectrum comes from the modes with low eigenvalues of DD'^ and 
D^D. In the strong magnetic field limit, we see from Figure [5] that the zero modes and the near-zero-modes all have 
spin up, (T3 = -|-1, as expected. The dominant contribution to the electric and magnetic moments are therefore: 

— / 771"^ \ / 771^ \ 

m(V^Ei2V) « tr2x2 --^T^ + tr2x2 -^^^^ (79) 



m \ I m 



™(^E34^) « -tr,x2 ) + tr2x2 ^,^^2^^ ) m 

For the magnetic dipole moment, the main contribution comes from the zero modes, so we simply count the degen- 
eracies in the various sectors: 



m(i/'Ei2V') 



B-F \ /'F_\ ( B + F \ / 

which is linear in the magnetic field B. For the electric dipole moment, the near-zero-modes cancel, leaving just the 
zero mode contribution: 

= ^ (82) 



which is independent of B, and negligible compared to BF, for B ^ F. [Note that (82) does not imply that there is 
a residual electric dipole moment when B vanishes, because (82) applies only in the B ^ F limit.] Thus, we see that 
the zero modes and near-zero-modes imply that 

(T^EiaV) « B , (ViEi2^) » (i/iE34V) (83) 



15 



This is in agreement with the lattice results of |26) . 

If we now consider the fluctuations in the electric dipole moment, we find a dependence on B, because 

(■(/; 1:341/' ^^34-4^) = tr ( -T-^^ — E34 -7-^ E34 

\J> + m ip + m 

_ / to2 n nt 1 

/ 1 ^ 1 

\W^WW ' ^ {m^ + DWf'im^ + DWf 

where in the last step we have used the fact that the dominant contribution comes from zero modes and near-zero- 



modes, all of which have (73 = +1. Thus, comparing with (81) we see that the fluctuation is linear in B 



(VS34V'V'S34^) « ( , ) B (85) 



27r2m2i4 



again in agreement with the lattice results of 



V. SMALL INSTANTON LIMIT 

In the opposite limit of a weak magnetic field the radius of the Landau orbit is much larger than the instanton size: 
1/y/B 3> p. It appears that a quantitative analysis is more difficult in this case, because all Landau levels contribute. 
However, we can still outline the qualitative picture. The effects induced by the presence of a small instanton on 
quark dynamics can be described in terms of the effective lagrangian introduced by 't Hooft (6] : 

C{x) = K e'" det [-^r{x)^Pl{x)] + h.c, (86) 

where k is a constant that contains exp{—8TT'^/g^), 6 is the 6'-angle of QCD (that we will assume be equal to zero), 
and the subscripts L and R refer to the left- and right-handed quark helicities. The flavor determinant breaks the 



UL{Nf) X Uii{Nf) symmetry of QCD down to SULiNf) x SUii{Nf). The (anti)instanton vertex as described by ([86]) 
absorbs Nf left-handed fermions and creates the same number of right-handed ones, or vice versa. 

Let us now turn on an external magnetic field directed along 2:3, and consider the electric dipole moment given by 



( 75 ) . The Wigner-Eckart theorem tells us that the only possible orientation of the dipole moment is along the 0:3 



axis. Using the chiral representation of the quark spinors as in (86), and using the gamma matrices (76), we can write 
down the electric dipole moment as 

O'f = -ipLOSipL + ^pRcrsipR ■ (87) 

This expression makes it clear that in the absence of an asymmetry between the left- and right-handed fermions the 
electric dipole moment should vanish identically, as required by P and CP invariances. However, the (anti)instanton 



transition caused by the interaction (86) can create a local asymmetry between the left- and right-handed fermions, 
and thereby induce a non-zero local electric dipole moment. 

It is interesting to discuss further the physical origin of this effect. In the absence of an instanton, in a magnetic 
field and with a fixed spin projection (73 = -1-1, there exist an equal number of left- and right-handed zero modes. 



as we discussed above in Section HA As a result, they cancel each other in (87). However, the instanton induces 
couplings between the fermion's spin, isospin (that belongs to the SU{2) subgroup of the color group SU{Nc)) and 
the (four-dimensional) orbital angular momentum, so that the operator contains the isospin-orbit term ^ T ■ L 
[B]. Because of this, the angular momentum of the fcrmion ceases to be a conserved quantum number ~ only a 
combination J = L + S + T \s conserved. The presence of the isospin-orbit term leads to the mixing of the lowest 
Landau level of left chirality with a radial excitation of right chirality that has the opposite parity. The change of the 
orbital angular momentum is compensated by a rotation in the color SU (2) sub-space. Because of this, the system 
acquires an electric dipole moment erf' ~ k signaling a local violation of parity invariance. The emergence of a local 
quark electric dipole moment in an external magnetic field has been observed on the lattice in |26j . where it was also 
observed that the electric dipole moment is strongly correlated with the local chiral density signaling the presence of 



the instanton, or some other topological object. Of course, at = the action (86) does not generate an asymmetry 



between the instantons and anti-instantons, so there is no global violation of parity. 
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VI. CONCLUSIONS 



While the physics of fermions in magnetic fields and instantons, separately, is well known and well understood, we 
have shown here that the combination of both background fields leads to a surprisingly intricate and rich structure 
in the Dirac spectrum. The inherent asymmetry when both instanton and magnetic field are present can lead to the 
development of an electric dipole moment. Physically, it can be understood as the outcome of two competing effects: 
the spin projection produced by a magnetic field and the chirality projection produced by an instanton field. We have 
illustrated this in detail both in the strong magnetic field limit in which the instanton scale is large compared to the 
magnetic length, and also in the opposite limit of small instantons. We have used the language of a four dimensional 
torus in Euclidean space, motivated in part by a desire to connect analytic results with recent lattice studies, which 
have shown a wide variety of interesting effects arising from the coupling of QCD to electromagnetic fields ;25.-.28, 37|. 
Corrections to the large instanton limit case could be constructed using the derivative expansion, in the natural 
Fock-Schwinger gauge [3Hj, x^^A^ = 0, which can be chosen for the combined instanton-magnetic field background, 
and which is an efficient means for computing induced currents, expectation values and correlators [39 41 . While 
in a constant self-dual field there is no preferred position at which the modes are localized, in the next order of the 
Fock-Schwinger gauge expansion we expect that these zero modes and near-zero-modes will localize on the instantons, 
as in the phenomenon of dynamical localization in the quantum Hall effect |42| . 

We thank T. Blum, M. Polikarpov and E. Shuryak for helpful discussions. This work was supported by the 
US Department of Energy under grants DE-FG02-92ER40716 (GB and GD) and DE-AC02-98CH10886, DE-FG- 
88ER41723 (GB and DK). 
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